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Abstract 

We study phase synchronization for a ratchet system. We consider the determin- 
istic dynamics of a particle in a tilted ratchet potential with an external periodic 
forcing, in the overdamped case. The ratchet potential has to be tilted in order 
to obtain a rotator or self-sustained nonlinear oscillator in the absence of external 
periodic forcing. This oscillator has an intrinsic frequency that can be entrained 
with the frequency of the external driving. We introduced a linear phase through 
a set of discrete time events and the associated average frequency, and show that 
this frequency can be synchronized with the frequency of the external driving. In 
this way, we can properly characterize the phenomenon of synchronization through 
Arnold tongues, which represent regions of synchronization in parameter space, and 
discuss their implications for transport in ratchets. 

Key words: Synchronization; Ratchets; Brownian Motors; Classical Transport 



PACS: 05.45.Xt; 05.40. Jc; 05.45.-a; 05.60.Cd 



1 Introduction 



The phenomenon of synchronization is widespread in Nature. We witness its 
manifestations in many different places and contexts. Synchronization is essen- 
tially a nonlinear phenomenon and is very common in many complex systems, 
not only in the physical sciences, but in the life sciences as well [1,2,3,4]. In 
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particTilar. the case of phase synchronization estabhshes a common formal- 
ism to treat both nonhnear periodic osciUators, as well as chaotic and noisy 
oscillators [1,5]. 

In a different context, there has been an increasing interest during recent 
years in the study of transport phenomena of nonlinear systems that can ex- 
tract usable work from unbiased non-equilibrium fluctuations. These, so called 
Brownian motors (or thermal ratchets) can be modeled by a Brownian parti- 
cle undergoing a random walk in a periodic asymmetric potential, and being 
acted upon by an external time-dependent force of zero average. The recent 
burst of work is motivated by both, (i) the challenge to model unidirectional 
transport of molecular motors within the biological realm and, (ii) the po- 
tential for novel technological applications that enables an efficient scheme to 
shuttle, separate and pump particles on the micro- and even nanometer scale 
[6,7,8,9,10,11]. 

Although the vast majority of the literature in this field considers the pres- 
ence of noise, there have been attempts to model the transport properties 
of classical deterministic ratchets as well [12,13,14,15,16,17,18,19]. In this pa- 
per we will be dealing with a deterministic tilted ratchet in the overdamped 
regime that acts as a rotator or self-sustained oscillator with a characteristic 
frequency, even in the absence of an external periodic forcing. The dynamics 
can be represented by a particle in a washboard potential that has been stud- 
ied in many different contexts, like phase dynamics in synchronization [1,5], 
pendulum dynamics [20], rotators [1], superionic conductors [21], optical po- 
tentials [22], excitable systems [23], diffusion on surfaces [24,25,26,27,28,29], 
charge density waves [30] and Josephson junctions dynamics [1,21,31]. When 
the washboard potential is periodically driven it exhibits a great variety of 
nonlinear phenomena including phase locking, hysteresis [32] and chaos [33]. 

Here we will study the synchronization properties of an overdamped particle 
moving on a tilted ratchet potential that is rocked by a periodic external force. 
Throughout this paper we will consider a constant force above the critical 
value, in such a way that the particle slides down the washboard potential, even 
though it is in the overdamped regime. When the periodic forcing is absent, the 
particle experiences only a fixed washboard potential, and moves through each 
period of the ratchet in a given constant time that defines the period tq of this 
rotator. The asociated frequency ujq = 271 /tq is its characteristic frequency. In 
this sense, this rotator is effectively acting as a self-sustained oscillator, with 
its own characteristic frequency. That is, if the rotator is driven by a constant 
force, it acquires the same features of a self-sustained oscillator, having a 
limit cycle in phase space. Thus, forced rotators are similar to self-sustained 
oscillators and can be synchronized by a periodic external force. We will drive 
this rotator with an external periodic force of period 0;^- In this way we can 
define properly the synchronization of the rotator and the external forcing. 
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The current or average velocity, which is the important quantifier for this 
system, displays steps as a function of a control parameter. This last result 
has been found previously by other authors that have studied overdamped 
deterministic ratchets [12,13,14,34]. On the experimental side, these so called 
Shapiro steps have been found recently for deterministic Josephson vortex 
ratchets and three- junction SQUID rocking ratchets [35,36]. 

More recent studies consider the problem of synchronization of deterministic 
ratchets, but they deal with complete synchronization between two coupled 
ratchets [37,38,39,40], and with anticipated synchronization between two uni- 
directional coupled ratchets with time delay [41]. In this work, instead, we are 
dealing with phase synchronization through a linear phase, properly defined 
through a set of discrete time events. 

As a way of characterizing the synchronization phenomenon, we will calculate, 
for the first time, the so called Arnold tongues for the tilted deterministic 
ratchet. For a description of Arnold tongues in circle maps and pendulum 
dynamics see [20,30,42]. Arnold tongues are regions of synchronization in a 
parameter space. Here we will calculate these regions in a two-dimensional 
parameter space defined by the ratio ujd/ojq and the amplitude of the driving 
periodic force Fj^. The tips of the tongues are located on rational values of 
the ratio LOj^jiOQ = p/q, where p and q are integer numbers. Each tongue is 
therefore labeled by a rational p/q whose inverse is precisely the value of the 
current in the driven washboard potential; the widths of these Arnold tongues 
correspond to the size of the steps of the current as a function of Fq. 



2 Tilted ratchets as nonlinear rotators 



To start out, let us consider now the one-dimensional problem of a particle 
driven by a periodic time-dependent external force in an asymmetric periodic 
ratchet potential. Here, we do not take into account any sort of noise, meaning 
that the dynamics is deterministic. Two additional forces act on the particle: 
a dissipative force proportional to velocity, and an external constant force. We 
thus deal with a rocked deterministic tilted ratchet [12, 14] in the overdamped 
limit that obeys the following equation of motion: 

dV(x) 

m-fx + ^ F + Fd cos{ujDt) , (1) 



where m is the mass of the particle, 7 is the friction coefficient, V{x) is the 
asymmetric periodic ratchet potential, F is a constant force, Fjy and repre- 
sent the amplitude and the frequency of the external driving force, respectively. 
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The ratchet potential is given by 

V{x) - Vo 



^ . 27r(a; — Xo) 1 . Anix — xq) 

6 — sm sm — ^— 

L A L 



(2) 



where L is the periodicity of the potential, Vq is the amplitude, and C is an 
arbitrary constant. The potential is shifted by an amount Xq in order that the 
minimum of the potential is located at the origin [16]. 

Let us define the following dimensionless units: x' — x/L, x'q = Xq/L, t' = 'jt, 
co'jj = a;B/7, F' = F/mL-f'^, F^ = Fo/mL-f'^, V = V/mL'^-i'^ and = 
Vo/mL^7^. Thus, we are using the periodicity of the potential L as the natural 
length scale and the inverse of the friction coefficient 7 defines the natural time 
scale. With these two quantities, the natural force is given by mL^'^ and the 
associated energy by mL^7^. 

The dimensionless equation of motion, after renaming the variables again with- 
out the primes, becomes 



x + 



dV{x) 
dx 



F + Fn cos(t<Jz)t), 



(3) 



where the dimensionless potential can be written as 



V{x) = Vo 



C — sin 2n{x — xq) — - sin A7r{x — xq) 



(4) 



and is depicted in the inset of Fig. 1. The constant C is such that V^(0) — 0, 
and is given by C = — (sin27ra;o + 0.25 sin 47ra;o). We choose, Xq ~ —0.19 



We can rewrite the equation of motion Eq. (3) as 

dU{x,t) 



x + 



dx 



0, 



(5) 



where U{x,t) — V{x) — [F + cos(a;£)t)]x. The important point to stress 
here is that we need to add a constant force F to the ratchet in order to tilt 
the ratchet potential and, in this way, obtain a rotator (that acts as a self- 
sustained oscillator), even without the external periodic forcing [1,44]. In this 
way, we can properly synchronize the characteristic frequency of the rotator 
with the driving frequency cuo- 

When = 0, we have a tilted ratchet that obeys the equation of motion: 
X + dV{x)/dx = F. The tilted (time- independent) washboard potential is, in 
this case, U{x) = V{x) — Fx, see Fig. 1. Thus, this ratchet becomes a rotator 
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Fig. 1. The tilted washboard potential indicating the dynamics that defines the 
discrete time events tk- The ratchet potential without tilt is illustrated in the inset. 

that has a characteristic frequency ujq. The associated period of the rotator 
To = 2tt/ujo can be obtained directly by integrating this equation of motion. 
However, there is another way to obtain uq for this tilted ratchet that relays 
on the introduction of a phase variable for this rotator. So, in what follows we 
will introduced this general concept that we will use in the rest of the paper. 

In order to define a phase variable we need first to obtain a discrete process 
from the continuous dynamics by introducing discrete time events. These dis- 
crete times can be defined as the times when the particle arrives at the discrete 
position Xk = ^k, which correspond to the minima of the ratchet potential 
without tilt. Here k = 0, 1,2, .... Remember that the period of the ratchet is 
one: V{x + 1) = V{x). This defines the set of times t^, where is a nonnegative 
integer. In Fig. 1 we show the washboard potential illustrating these discrete 
times. Once we obtain these set of markers, we can define an instantaneous 
linear phase for the rotator as [1,5] 

(t>{t) = 2tt /'^\ + 27r/c, h<t< 4+1 (6) 

This linear phase is valid in the indicated interval and defines a piecewise 
linear function of time that increases by 27r each time the particle crosses the 
dimensionless position Xk = ±k. 



5 



Given this phase, we can define the instantaneous frequency of the rotator as 
u!{t) — and the average frequency as 

T 

{00) = hm 1 fu;{t)dt = hm - 0(0)]. (7) 

1 — >oo 1 J 1 — >oo _Z 





Without loss of generality, we choose to — 0, and thus 0(0) — 0. The limit 
above can be written as 

(^) = lim ^ = 27r lim - (8) 

fc— >0O tf; fe— >C» tj~ 



This is the simplest way to calculate the average frequency; wc simply count 
the number of jumps (given by k) and divide by the total time span t^. 

The other quantity of importance is the average velocity or current in the tilted 
ratchet. In order to evaluate this current we have to calculate the number k 
of unit periods that the particle crosses to the right, denoted by N^, and the 
number of crossings to the left, given by N^. The total number of periods 
traversed on the ratchet is given by N'^ — + N^. The difference 

N, = - (9) 



indicates that during the time tk the particle has covered the distance x^ — Nj^. 
Therefore, the average velocity (current) is given by 

{v) = hm ^ (10) 



In the particular case when = 0, that is, when there are no jumps to the 
left, we have Nk = iVf = k. Thus, 

(v) = hm ^ = i^(^) (11) 

k^oo tf; ZTT 



In the simple case of a tilted ratchet without external forcing, Fd = 0, the 
average frequency defined above coincides with the natural frequency of the 
rotator, that is, (cu) — cuq. 

The above treatment is quite general and can be used in the case of a tilted 
ratchet with an inertial term, even though this inertial ratchet can display 
a chaotic dynamics [15,16,17,18,19], and also in the case of a tilted ratchet 
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Fig. 2. The average velocity as a function of the external force F. The dashed line 
indicates the case when the periodic driving is absent (Fq = 0) and the continuous 
line is the case when the periodic driving is acting on the particle with an amplitude 
Fd = 0.5. In both cases we used lod = 0.7. 

with noise. We have used the concept of a linear phase due to its broad-range 
applications to the cases of periodic and chaotic oscillators [1], chaotic rotators 
[44] , and oscillators in the presence of noise [5] . Additionally, the introduction 
of the discrete events that define the linear phase allow us to simplify the 
dynamics and have a more clear picture of the synchronization involved. 



3 Numerical results 

In this section we will solve numerically the equation of motion for the rocking 
tilted ratchet. We use the fourth-order Runge-Kutta algorithm to solve the 
differential equation Eq. (3). Once we obtain the full trajectory, we identify 
the set of discrete times tk when the particle crosses the positions Xk- With 
this marker events we calculate directly the average frequency using Eq. (8) 
and, after calculating the quantity Nk, we obtain the current, using Eq. (10). 
We will fix throughout the paper the amplitude of the ratchet potential as 
Vo = l/27r. With this value, the critical tilt to the right is = 0.75 and to 
the left F^^-1.5. 

In Fig. 2, we depict the average velocity, scaled with the driving frequency, 
27i{v)/u!d as a function of the tilt F. The dashed line shows the case without 
periodic driving (Fq = 0) and corresponds to the fixed washboard potential 
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Fig. 3. The scaled average velocity 2ti{v) /i^o as a function of the amplitude of the 
periodic driving Fj:,. (a) The dashed line indicates the case when the tilt is absent 
[F = 0) and the continuous line is the case when the tilt is = 1. (b), (c) and (d) 
correspond to successive magnifications of the current for the tilted ratchet, showing 
a self-similar structure of steps, typical of a devil's staircase. In the dashed line in 
(a), we used lod = 0.7, since in this case ujq is not defined. In the other cases we 
used ojd = O.Two- 

U{x) = V{x)—Fx. Notice that the current is zero until we arrive at the critical 
tilt = 0.75 to the right or to the critical tilt F^ = -1.5 to the left. This 
step of zero current is not centered around the origin, due to the asymmetry of 
the ratchet potential. For values greater than F^ we have a finite current that 
increases monotonically with F. Of course, for values less than F^ we obtain 
a negative average velocity that decreases for negative values of the tilt. When 
the periodic driving is present, this scaled current acquires a series of clearly 
defined steps for values of the current given by the ratio p/q, where p and 
q are integer numbers. In many cases, q = 1 and the average scaled current 
is an integer. In the context of Josephson junctions, these are the celebrated 
Shapiro steps [33]. 



Remember that in the case where all the jumps are to the right, that is, 
A^^ = 0, we show that 27i{v) = (u). Therefore, a rational value of 2'n{v) /ujd = 
p/q means that (a;) = {p/q)uiD for a whole range of values of the tilt. This 
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Fig. 4. Arnold tongues in the parameter space Fd against u)£,/ujq. Notice that the 
tongues are located precisely at the rational values p/q, for p and q integer numbers, 
of the ratio between the driving frequency and the natural frequency of the rotor. 
Here the tilt is F = 1 and the corresponding frequency is ojq ~ 3.41 . Each tongue 
is labeled by the inverse q/p that gives the value of the scaled average velocity 
2tt{v)/ujd in that region of the parameter space. 

phenomenon is called frequency locking. 

In Fig. 3, we show the scaled average velocity 2'k{v) /ojd as a function of the 
amplitude of the periodic driving Fr,. In (a), the dashed line depicts the current 
for the ratchet without tilt {F = 0) and coincide with previous calculations 
[12,13,14] showing a structure of steps of unit height. The solid line shows the 
current for a tilted ratchet with F = 1 that also has well defined steps for 
rational values. In (b), (c) and (d) we show successive magnifications of the 
current that clearly exhibits a self-similar structure of steps, typical of a devil's 
staircase [42,43]. This detailed structure has been reported before [13,14] for 
an overdamped ratchet without tilt, but here we obtained this fractal current 
with a devil's staircase also for the tilted ratchet. 

In Fig. 4, we depict the parameter space Fo against ujd/uq that shows regions 
of synchronization, called Arnold tongues, located at rational values p/q, where 
p and q are integer numbers. We choose a tilt F = 1, which corresponds to ujq ~ 
3.41. Notice that the tongues start, for small values of F^,, precisely at these 
rational values, as indicated in the figure. Each Arnold tongue corresponds to 
one particular rational p/q whose inverse gives the value of the scaled average 
velocity 2tt{v)/ujd in that region of the parameter space. This is depicted in 
a three-dimensional plot in Fig. 5. Therefore, the average velocity, properly 
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Fig. 5. Three-dimensional plot of the scaled average velocity 2'k{v) /lod as a function 
of Fd and ujd/^jJq- A projection of this 3D plot shows the Arnold tongues in the 
parameter space of Fig. 4. Here the tilt is F = 1 and ujq ~ 3.41. 

1-4 I 1 1 1 1 1 1 
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COq/COq 

Fig. 6. Average velocity, scaled now as 2tt{v)/loo, and the average frequency of the 
rotator, scaled as [uj) /ujq as a function of the ratio ud/^o - The straight line segments 
correspond to the steps observed due to phase synchronization in Fig. 4. The peaks 
in the current are located at the right borders of the Arnold tongues in Fig. 4. Here 
the tilt is F = 1, Wo - 3.41, and Fd = 0.5. 

scaled, acquires rational values g/p, which correspond to the phenomenon of 
phase synchronization in this forced tilted ratchet, acting as a rotator [44]. 
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Finally, in Fig. 6, we plot the average velocity, scaled now as 2Tr{v)/uJo, as 
a function of the ratio ujjj/ujq. Here cjq is a fixed value that correspond to 
the characteristic frequency of the tilted ratchet in the absence of periodic 
driving. In the same figure, we plot the average frequency of the rotator, scaled 
as {uj)/ujq, calculated using the discrete dynamics explained in the previous 
section. Instead of plateaus, we have now straight lines, since we are scaling 
the current with a fixed parameter c<Jo, instead of the running parameter ujd. 
In this case, all the jumps arc to the right direction and therefore 27r{v) = (a;), 
that is, the average velocity is proportional to the average frequency of the 
rotator. 

Notice that the peak values of the average velocity correspond to the borders 
of the Arnold tongues in Fig. 4. Thus, when we are crossing an Arnold tongue 
(synchronization region) the current increases linearly with ujd/oJq] at the right 
border of the tongue the current is maximal and outside the tongue starts to 
decrease until we arrive at the next tongue to increase linearly again, and so 
forth. At the right border of the Arnold tongue, labeled hj p/q = 1, the current 
has a maximum, followed by the second largest peak at the right border of 
the tongue with p/q — 2. The heights of the peaks in the current arise due 
to the combined effect of both the width of the tongues and the slope of the 
linear segments. The slopes of each of the segments correspond precisely to the 
inverse values q/p that label the tongues. Therefore, the peaks in the current 
are associated with the phenomenon of synchronization. 



4 Concluding remcirks 

In summary, we have analyzed the phenomenon of phase synchronization in 
tilted deterministic ratchets in the overdamped regime and with an external 
periodic forcing. The dynamics in this rocked washboard potential corresponds 
precisely with the dynamics of a rotator that, acting as a self-sustained oscil- 
lator, can be capable of being synchronized with the external periodic drive. 
We can clearly identify three frequencies for this system: the characteristic 
frequency of the rotator without driving, the driving frequency itself, and the 
average frequency of the rotator with driving. This average frequency is the 
derivative of a time-dependent phase, that can be obtained through a set of 
discrete time events and is a piecewise linear function of time between these 
markers. We calculated the average frequency and the average velocity as a 
function of the tilt and obtained the well-know Shapiro steps that characterize 
the phenomenon of frequency locking. We also exhibit the self-similar structure 
of steps in the current, typical of a devil's staircase. We obtained well-defined 
Arnold tongues in the 2D parameter space given by the amplitude and the 
frequency of the periodic forcing. Each Arnold tongue is labeled by a rational 
number p/q, where p and q are integer numbers, whose inverse gives precisely 
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the rational value of the average scaled velocity of the particle. Finally, we 
show that the local maxima in the average velocity correspond to the borders 
of these Arnold tongues and, in this way, we established a connection between 
optimal transport in ratchets and the phenomenon of phase synchronization. 
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